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Abstract: In this paper we study the sine-Gordon and the Liouville hierarchies in laboratory coordinates 
from a bi-Hamiltonian point of view. Besides the well-known local structure these hierarchies possess a 
second compatible non-local Poisson structure. 
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O ; 1 Introduction 
0\ 
O 

' In [1] Dubrovin and Zhang started a classification program of integrable evolutionary PDEs based on the 
I bi-Hamiltonian geometry. 

> 

X 



As a first step in the realization of this program they studied bi-Hamiltonian hierarchies of PDEs of 
the form 

oo 

ui = Ai{u)ui +Y,e''Pk{u,u^, ...) (1.1) 

k=l 

where Pk{u,Ux, ■■■) are differential polynomials. 
A central role in their approach is played by the bi-Hamiltonian structure of the dispersionless limit 
(e = 0) of the hierarchy (|l.lj) . Such structure consists of a pair of compatible local Poisson brackets of 
hydrodynamic type (see 0). 

Very important examples of integrable PDEs such as the KdV equation, the continuous limit of 
the Toda lattice and, after a suitable change of coordinates, the NLS equation belong to the class of 
hierarchies considered by Dubrovin and Zhang. 

An open problem is to understand how to extend the classification scheme of |H onto more wide 
class of integrable systems not admitting dispersionless limit. The first point to clarify is which kind of 
bi-Hamiltonian structure is involved in this case. 
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In this paper we condider two important examples of this class of systems: the sine-Gordon equation, 
one of the first PDE solved by inverse scattering transform (see ^ and [S]) and the Liouville equation: 



iott - uj XX = -sin LO (1.2) 
(ptt - 4>xx = -e"^ (1-3) 

(in laboratory coordinates). 

Starting from the results of a paper of Boiti,Leon and Pempinelli ([2]), we show how to construct 
their bi-Hamiltonian structures. 

In sine-Gordon case such structure consists of the following compatible pair of Poisson bivectors 



Pi 



-2i 



(1.4) 



and 



'Sl,^d-'SI,^ + + s)d + dii + .) - f Sl,^d-'Sf,^ - 




(1.5) 



with Sj-Q^ = ^ ( t ^ + iS^Q^ = —2u and s = sq cos uj , while in the Liouville case we obtain the pair 



Pi 




[1.6] 



and 



I fd- S[l^d-^S[l^ + (8^2 + e'^)d + d{8u^ + e<^) - Ad^ -S[l^d-^S[l^ + 32u2 + 2e<^ + 2 - 89^ 
16 1^ -s[l^d''S[l^ - 32n2 - 26-^ - 2 + 8d^ -S[l^d-'S[l^ + 169 

where = ^ + Au^ and S[l^ = -8u. 

Following the Lenard-Magri scheme (\i-2\) one can construct recursively the integrals of motion. In the 
sine-Gordon case these integrals usually appear in two distinct families corresponding to the asymptotic 
expansion of the reduced monodromy matrix T(A) around A = and X = oo { 7 ) .In the bi-Hamiltonian 
framework these two families are arranged in a single hierarchy, infinite in both directions. 

Our observation is that, in both cases, the bivectors H1.5|l and H1.7|l are weakly nonlocal according 
to the definition of Maltsev and Novikov (see [ri| ) . 
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This class of nonlocal Poisson bivectors appeared only recently in the mathematical literature as a 
generalization of certain nonlocal Poisson brackets of hydrodynamic type (see jH] , J3] and ^1] ) . 

Nevertheless the list of (1 + l)-integrable systems admitting such kind of Poisson structures contains, 
besides the sine-Gordon and Liouville equations, many other well-known integrable systems as KdV 
(considering higher Poisson structures) and NLS (see [llj'). 

To the best of our knowledge the bi-Hamiltonian structure of the Liouville hierarchy was unknown, 
while the bi-Hamiltonian structure of the sine-Gordon hierarchy appeared the first time in the paper [H] 
(see remark 5 in section 3). 

2 Preliminary definitions 

Let L be the space of maps from the real line to some n-dimensional manifold M endowed with local 
coordinates ti^, and let A be the space of differential polynomials in n*''', that is 

f^A^f = j;/n,si,.,.„,s„(^x)n(*i'^i)...u(*™'^-), (2.8) 

where 

n(*'^) = < 

The role of the functions on C is played by the local functionals 



" + 00 

/ = / f{u{x),Ux,Uxx, ■■■)dx 

-oo 



where f ^ A. 



Definition 1 A (non local) multivector X is a formal infinite sum of the type 
X = X'^^'''''-'"""^{xi,...,Xk;u{xi),...,u{xk),u^{xi),...) ^ . ^ . , A ... A 



where the coefficients satisfy the skew-symmetry condition with respect to simultaneous permutations 

^p: Xp <— > Iq, Sq, Xq 

Definition 2 A k-vector is called translation invariant if 

X^^^'-^'-'^'""'{xi,...,Xk;u{xi),...,u{xk),...) = a^i...a^^X*i'-'^'=(xi,...,Xfc;n(xi),...,u(xfc),...) 
where X*^'''''*'''(...) means X^'^'^'''''^'"^{...) and for any t 

X'^'-'"'{xi+t,...,Xk+t;u{xi),...,u{xk),...) = X'^'-'^^ixi, ...,Xk;u{xi), ...,u{xk), ...) 
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It follows from this definition that a translation invariant multi- vector field is completely characterized 
by the "components" 

Xii-ik .- X{xi, ...,Xk;u{xi), ...,u{xk), ...) 
Definition 3 A bivector P is called Poisson bivector if and only if 

[P^ P]schouten — (2-9) 

We recall that the Schouten bracket of two translation invariant bivectors is given by the formula (see 



0) 



rp (^]ijk _ 1 ( dPxy r^sf^lk I dQiy If. dPxy r^Sf^lk , ^Qjy 

apki riO^^ fipki FtO^^ flP^^ 



A Poisson bivector defines a Poisson bracket. Indeed we have the following 
Definition 4 The Poisson bracket of two local functionals Ii ,I2 is given by the formula 

Definition 5 A pair of bivetcors (Pi,P2) defines a bi-Hamiltonian structure if and only if 

[Pi, Pi] Schouten = [Pi, P2] Schouten = [P2, P2] Schouten = 0. (2-12) 

Definition 6 A Poisson bivector P of components Pxy (and the corresponding Poisson bracket) is called 
local if it can be written as: 

P% = Y.B'i{u,Ux..-)5^'\x-y) (2.13) 
k 

By substituting (^T^ in we get 

^ ^' J 6u^{x) 6u^{x) ^ ^ 

where 

^*'=E^fc'("'^---)^' (2-15) 



Definition 7 A Poisson bivector P of components Pxy (and the corresponding Poisson bracket) is called 
weakly non-local (see ',11]) if it can be written as: 

P'Jy = Y,B'^{u,Ux,....)5'^^\x-y) + Y,euSl^^^^^ (2.16) 

where v{x — y) is the step function (= ^ ifx — y is positive and = — ^ ifx — y is negative). The constants 
Gkl = &lk define a quadratic form in the space generated by the flows S^^j^-^{u,Ux, •••)• 

By substituting (|TT^ in (|TTT|) we get 

{/i,M:=p-^P-^dx (2.17) 

where 

P'' = Yl ••••)^'' + E ■■■)d''si^^{u, Ux, ...) (2.18) 

k k,l 

and, by definition 

^"'(•):=o / -/ ]dx{.) (2.19) 



2 \ / / / 

Later on we will call local and weakly non local bivector also the corresponding differential and 
integro-differential operator. 

Similarly it is possible to define the class of weakly nonlocal symplectic structures (see [TT]'). 

Definition 8 Let P be an invertible Poisson bivector. Then P^^ is a weakly nonlocal symplectic structure 
if and only if it can be written in the following way: 

= + Z]^'^Afc)(^'^^'---)^"^'3i)(?^,?^x,..^ (2.20) 

3 The sine- Gordon hierarchy 

3.1 Useful results 

• The starting point of the paper P] is the hierarchy of equations represented by the equation of zero 
curvature: 

Ut-Vx + [U, V]=0 (3.21) 



where (see P] for details) 



U = -iXa3+ucri + {^){sa3+iva2) (3.22) 
A 

n p 

V = Y^ VjX"-^ + ^ VFj A^-f-^ (3.23) 
i=o i=o 
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aj are the Pauli matrices and the functions u{x), v{x) and s{x) have the following behavior for large 

x\: 



/im|^.|_ 


^oou{x) = 




^oov{x) = 


limM_ 


^ooS{x) = So 



The sine-Gordon hierarchy corresponds to the reduction of this hierarchy on the surface 

s2 - ^2 ^ si (3.24) 



The equations of this hierarchy can be written as 

r(+) k 



'i-'E^r'.^^i+'E^r'^^i ,3.25) 



where 



j=0 j=0 



\ov S OS / 



and 



(3.26) 



The first "Hamiltonians" are: 



oo 

oo 



J = ^ M 3.27 

\-2six) J 



dx{s-so + -u^] (3.28) 



oo 



h[ ^ = -i I dx { s - So + ^W^ + ^ + 4.u) ] (3.29) 



The sine-Gordon equation corresponds to the choice m = k = 0, fi^^^ = '* ~ ^' t^y~' ~ /^j = 
for J different from 0. In fact by straightforward calculation one gets: 

(3.30) 



.Vtj \-Vx-AuSj 



On the surface (13.241) v and s can be written as functions of 



ijj: 



V = ISO sin w (3.31) 
s = So cos a; (3.32) 
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By inserting (|m-{2j) and (|H.H2|) in the system one obtains sine-Gordon equation 

uJtt - ^xx = -16so sin w 



(3.33) 



The recursion operators for this hierarchy are given by the formulae 

N -- 



{2v + u^)d-^u + + s - f {2v + ^x,.)5"'f + T - f 



-2sud + \sd 



(3.34) 



r,2 



(3.35) 



si \\s {-d + wd-^w) swd'^ (f + \wd\) - \sd^\ + s ) 

where w = ^ + 2u. 

• The hierarchy is local, that is all the equations are differential. 
3.2 Prom the recursion operator to the bi-Hamiltonian structure 

In this section we show how to construct the sine-Gordon bi-Hamiltonian structure starting from the pair 
{J,NJ), where NJ = 



-{2v + u^)d-^{2v + n^) + + s)d + d{W + s 



2' 



91 
4 



92 



Y 2usd-^{2v + lix) - 2s{u^ + s)+ ^ 

We begin with the following 
Remark 1 The system 



2{2v + Ux)d~^us + 2s{u'^ + Sj ^ 
—4sud~^us + sds 



(3.36) 



P21 P22, 



m 

5u 

SH I V SH 
^ 5v s 5s / 



is equivalent, on the cylinder \3.24-^ , to the Hamiltonian syst 



em 



(3.37) 



In fact the third equation 







P12 


A2f \ 


ou 


Vt 


= P21 


P22 


p V 
P22-, 


SH 

Sv 


\st) 




V p 


s 22 g y 


\Ss/ 



V „ 6H V „ 6H w V 6H v 

St = -P21-, h --P22^ + --P22 — r~ = 

s ou s ov S S OS s 



is compatible with the constraint 



(3.38) 



(3.39) 
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Then we have to substitute the 2x2 matrices J and N J with the corresponding 3x3 matrices J = 

^ Jii Ji2 Ji2^ ^ 
J21 J22 J22f (3.40) 

\-^'J21 ^J22 7 -'22 7/ 

and NJ = 



( {NJ)u {NJ)i2 {NJWj ^ 

{NJ)21 {NJ)22 (iVJ)22f (3.41) 
\^ANJ)2l f(iVJ)22 f(iVJ)22f/ 

The idea is that the sine-Gordon bi-Hamiltonian structure is the restriction of the pair (J, NJ) on 
the surface 1)3 .241) . By replacing the coordinates {u,v,s) by the coordinates {u,u;,so) where 



s = sq cos uj 



(3.42) 
(3.43) 



and taking into account that, after a change of coordinates 



v}' (u^) a bivector P^^ transforms as 



we obtain 



and 



(l 





AS ■ V 

so so 



J21 







kl 














J12 






A 





\ 


(d 


-2i 


o\ 


-^22 


T ^ 







■^0 


y_ 

So 


= 2i 








V 7 

7 -'22 


V J V 

7 "^22 7/ 






■^0 


so / 







0/ 



(3.44) 



/l 


\ 


1 {NJ)n 


{NJ)l2 


(iVJ)l2f 


\ 


A 





\ 





AS ■ V 

^0 *0 


{NJ)21 


{NJ)22 


(iVJ)22f 







*o 


V 
so 




y_ _s_ 

So So / 


\^{NJ)21 


%NJ)22 


%NJ)22\ 


) 




4 


so / 



'-{2v + u^)d-^{2v + u^) + (^ + s)d + + s) 
-2iud-^{2v + u^) + 2i(u2 + s) - ^ 

V 



91 
4 



-2i(2u + ti^)a"in-2i(n2 + s) + ^ o\ 
4ud~\ -d 

Oy 



^From these formulae it follows immediately that the restriction on the surface (|3.24l) is given by the 
expressions (|1.4() and ()1.5|) . 

Then we have the following theorem. 
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Theorem 1 The pair (J, NJ) and its restriction on the cylinder \3.24\) are bi-Hamiltonian structures. 



We omit the proof of this theorem consisting of a very long but straightforward calculation based on 
the formula (|TTn|) . 

The most difficult part of this kind of computation is to check whether certain expressions containing 
products of distributions vanish. The trick is to use some identities between distributions like 

- y) = V ( ' ) f'^''\x)6('-''\x - y) 
6{x — y)6{x — z) = 6{y — x)5{y — z) = 5{z — x)5{z — y) 

v{y — z)5{z — x) = dy^ {6{z — x)6{z — y)) = dy^ {5{x — y)6{x — z)) = —i>{x — y)5{x — z) 
v{x - y)6{y - z) = d^^ {6{y - x)6{y - z)) = d^^ {6{z - x)6{z - y)) = -v{z - x)5{z - y) 
v{z — x)5{x — y) = d^^ {5{x — y)S{x — z)) = d^^ {5{y — x)5{y — z)) = —5{y — x)v{y — z) 

^From the theorem it follows directly 
Corollary 1 M.5\) is a weakly non-local Poisson bivector. 

Remark 2 The two hierarchies {In+i} '■= {~^-f^2n+i} '^"-^ i^-n} '■= {^^2n+i} '^^'^ "glued" together. 
In fact we have: 

P^dlo = —iP26H+ = -iPi6H- = Pi6Ii (3.45) 
So 
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In other words we have the Lenard-Magri chain: 



1-2 



h 



Pi 
/ 



Pi 
\ 



Pi 



Pi 



Pi 
/ 



Pi 
\ 



Pi 
/ 



Pi 
\ 



Pi 
/ 



Pi 
\ 



P2SI0 = PiSh 



P2SI1 = P1SI2 



An immediate consequence of the last remark is the following 

Theorem 2 The sine-Gordon equation belongs to the hierarchy {Hj^ = sqI^ + Ik+i}kez- 
where 

,2 



HisG = j [u^ + ^u^x - + 2so cos a; ) dx (3.46) 

H2SG = 



/■ 1 A 2 1- \ ^ ( 



' ° cos(2w) - — ^ COS uu^ - + ■ 



V 4 8 32 ' 128 

Sq 2 ^ -^2^'^^22^'^ "^417/^ i-r\ 

+—U COSCJ + lUSo COSUUx - -^ULOxxx + o«a; " TtULO^ - —U UJ^ + -U U}x + -u dx (3.47) 

z o o lb ib 4 o y 

The momentum integrals are given by the following expressions: 
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PiSG = -\j (4iwa;^ - ^l) 



dx 



(3.48) 



P2SG = 



1 /1 2 1 



1 fs, 



— ( 2" + 2''''^^ - Y + *o cosw ) - -2 { -fcos{2uj) -—cosiOiui-^ + j^ 



3so 



Sq \ 4 



So 2 



32 128 



1 



cosw + iusocosujujx - -^uiOxxx + g% ~ iq^^ ~ ig" '^^^ + -^u ^^x + -^u ] dx (3.49) 



Moreover the Poisson bivectors Pi and P2 are invertible. More precisely: 

'0 



and P^ 



— 



i _^ 

v2 4/ 



(3.50) 



/ Tid-^Ti + ld ri5-ir2-z(f cosa;-f -5) \ 
\T2d-^Ti + z (f coso; - ^ - T2a-iT2 + - f cosa;J 9 + ^ ^ ^ / 



(3.51) 



= - (za;^ + 2u) 



zsosino; - -^'^xx - -^x 



Proof By adding the Lenard-Magri chain for the Hamiltonians {/fejfegz multiphed by sq to the same 
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chain shifted by one we obtain 



Ho 



Pi 



Pi 



Pi 
/ 



Pi 
\ 



Pi 
/ 



Pi 
\ 



Pi 
/ 



Pi 
\ 



Pi 
/ 



Pi 
\ 



P2dHo = P16H1 



P2SHI = PISH2 



In order to find the recursion operators of this hierarchy we need the following 

Remark 3 Let N and M = be the 2x2 recursion operator introduced above. If M and N satisfy 
the equations 





Jl2 


J21 


J22 



Jll Jl2 
I "^21 "^22 y 



(3.52) 



and 



Jh 


Jl2\ 




J12 


^21 


J 22/ 




J 22 



(3.53) 
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then the 3x3 matrices 



and 



satisfy the equations: 



and 



where 



and 



M 



N 



1 Mil Mi2 \ 

M21 M22 

[^M2l fM22f/ 

^ Nil N12 \ 

N21 N22 

VfiV21 fiV22f/ 





MJ2 


= Jl 








= J2 








'^12 


•^12 s 


Jl = 


71 


•^22 


jl V 

•^227 




L /l 


1^22 


V jl V 

7'^22I/ 




-'11 


72 


7-2 H \ 
-^12 s 


J2 = 


"^21 


72 
•^22 


t2 1) 
•^22 7 




\ 5-^21 


V j2 
1^22 


V j2 V 
7 "^22 7/ 



(3.54) 



(3.55) 



(3.56) 



(3.57) 



(3.58) 



(3.59) 



(in our case Ji = ( [X^ and J2 = ( [JT^ j 

Since Ji and J2 are not invertible in the space {u, v, s) we can try to express M and N in terms of them 
only after the restriction on the surface (|3.24jl . 



Taking into account that after a change of coordinates 



M* (u^) M transforms according to the 



rule 



-Mr 



it is easy to see that in the coordinates (u, 01,50) M becomes: 

/ s — vd~^ 



(3.60) 



w 



|s9 - 2ivd-^ {v + Iwd) 
^[-d + wd^^w) s - \d'^ + wd~^ [v + \wd) 
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The equations (|H.5()j) in these coordinates becomes 





s — vd 


- 2ivd-^ 


{v + \wd) 




/ pii 


pl2 


o\ 




pl2 


o\ 


i 
2 


[—d + wd~^w) 


s - \d'^ + wd~ 


1 {v + Iwd) 


* 


p21 


p22 
^2 







p22 





V 










V 


\ 





0/ 


V 





0/ 



^From these equations it fohows immediately that 



M' 



s — vd 



^sd - 2ivd-^ {v + Iwd) 
si {-d + wd'^w) s - + wd~^ {v + ^wd)^ 



satisfies the equation 



M'P2 = Pi 



(3.61) 



and therefore it is a recursion operator for the sine-Gordon hierarchy. 

It is easy to check that the bivector p.5U() is the inverse of Pi and that the bivector (|1T3T|) = Pf ^M' 
is the inverse of P2. 
Finally from 



P16P1SG = P2SP2SG 
PiSHisG = P2SH2SG 



(3.62) 
(3.63) 



it follows immediately 



SP2SG = Pi^PlSPiSG 
6H2SG = P^^Pl^HisG 



(3.64) 
(3.65) 



By straightforward calculation we get 

5H2SG -1 . 3 2 3. 2 1 3 1. 1. 

J = —Uxx + ISq COSLOUJx - -ULO^ - -lU LO^ + -U + -ISQiOx + USq COS W + USq - -lUJxxx + 

ou 4 o 4 z z o 

1 , 



16 
3 



-isoCOSOJUx — -Sq sin2u; + -iuuJxOJxx + --:UUxOJx + -^u^^xx — Sq sin a; + 
2 8 4 8 



+ -sqcosujuJxx + ^sqlOxx " ]^isQUx + ^iUxxx " ^Sosinwu;^ + ^iuxujl - ^JluJxx + ^iu^Ux + 



1 



1 



--r^oJxxx - 7;U sosma; 
lb 2 



that is the differential of the Hamiltonian (|3.47j) . Similarly we get (|3.49|) . 

Corollary 2 The operator i'j.51\) defines a weakly nonlocal symplectic structure. 
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Remark 4 In order to obtain the usual symplectic structure of the sine-Gordon hierarchy 



Gi 



( 



-1 







1 



(3.66) 



it is sufficient to make the change of variables q = to, p = Aiu — qx: in the new coordinates the Poisson 



Fuchssteiner and Oevel (^). 

Remark 5 The geometrical structures (recursion operator, bi-Hamiltonian structure) of the sine- Gordon 



Starting from these structures it is not obvious how to get their counterpart in laboratory coordinates. 

A possible way was suggested by Fuchssteiner and Oevel in (^91 ) and it is based on the properties of a 
particular class of time- dependent vector fields called master symmetries. In the case of integrable PDEs 
admitting a hamiltonian formulation with respect to some Poisson bivector P a master symmetry r has 
the property that the Lie derivative LiCrP defines a new Poisson bivector compatible with P. 

The main idea of the paper jO^ is that the master symmetry of the sine- Gordon equation in light- cone 
coordinates can be modified in such a way that the new vector field coincides with it on the submanifold 
of the solutions and gives rise to an evolution symmetry with respect to the laboratory coordinates. 

Such evolution symmetry is the master symmetry of the sine- Gordon hierarchy in laboratory coordi- 
nates. Therefore the Lie derivative of the symplectic structure it^. 6'6|) along it defines a second compatible 
Poisson bivector @2 CL^d the operator <I> = 020r^ hereditary symmetry, in other words for any 

vector fields A and B it satisfies the equation: 



The proof of the hereditaryness of the operator $ is based on some computer algebra methods introduced 
in the paper .10^ . 

4 The Liouville hierarchy 
4.1 Useful results of |j2j 

• The Liouville hierarchy corresponds to the reduction of the hierarchy of equations 1)3.21(1 on the plane 
in the space (n, v, s) given by the equation s = v. 

• The difference with respect to the sine-Gordon case is that the recursion operator (|3.34|) is no longer 
invertible and as a consequence only the local Hamiltonians {-ff^jjgAr and survive. 



bivector \1.4\ becomes the usual one while the second one coincide with the Poisson bivector found by 



equation are well known in light- cone 




^\A, B] + [^A, ^B] = $ {[^A, B] + [A, ^>S]) 



(3.67) 
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4.2 The bi-Hamiltonian structure of the Liouville hierarchy 

We have seen that the bi-Hamiltonian structure (|1.4| 1.5j) of the sine-Gordon hierarchy is simply the 
restriction of the pair {J,NJ) on the surface — = Sq. Similarly it is possible to obtain the bi- 
Hamiltonian structure of the Liouville hierarchy as restriction of the pair (J, jqJ + NJ) on the plane 
s = V = The only difference in the Liouville case is that, in order to have a second local hamiltonian, 
we have to choose jqJ + NJ instead of NJ. 

Acknowledgments I would like to acknowledge Boris Dubrovin for useful discussions and for the 
reference [5], the starting point of the present work. I thank also Maciej Blaszak for the reference jH]. 
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